Fast magnetic reconnection in the plasmoid-dominated regime 
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A conceptual model of resistive magnetic reconnection via a stochastic plasmoid chain is proposed. 
The global reconnection rate is shown to be independent of the Lundquist number. The distribution 
of fluxes in the plasmoids is shown to be an inverse square law. It is argued that there is a finite 
probability of emergence of abnormally large plasmoids, which can disrupt the chain (and may be 
responsible for observable large abrupt events in solar flares and sawtooth crashes). A criterion for 
the transition from magnetohydrodynamic to collisionless regime is provided. 
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Introduction. Magnetic reconnection is the process of 
topological rearrangement of magnetic field, resulting in 
a conversion of magnetic energy into various forms of 
plasma energy [![. It is believed to cause solar flares 
and has been studied in tokamaks 0, dedicated labora- 
tory experiments Q and measured in situ in the Earth's 
magnetosphere The basic conceptual underpinnings 
of the modern understanding of resistive reconnection 
can be summarised in three points: (i) generic X-point 
configurations are unstable and collapse into current lay- 
ers [f| 111 ; (ii) the structure of resistive current layers is 
well described by the Sweet-Parker (SP) model Q: if 
Bo is the upstream magnetic field, Va = Bo/y/Arrp is 
the Alfven speed (p the plasma density), L the length of 
the layer, n the magnetic diffusivity, and S = VaL/t] the 
Lundquist number, then the layer thickness is S ~ Lj \[S, 
the outflow velocity is Va, and the reconnection rate is 
cE ~ VaBq/VS — "slow" because it depends on S, 
which is very large in most natural systems; (iii) when 
S exceeds a critical value S c ~ 10 4 , the SP layers are 
linearly unstable [§] and break up into secondary islands, 
or plasmoids This fact has emerged as a defining 
feature of numerical simulations of reconnection as they 
have broken through the S c barrier It seems 

that high-,!? reconnection generically occurs via a chain of 
plasmoids, born, growing, coalescing, and being ejected 
in a stochastic fashion 17 Jl8 |. Importantly, recent nu- 
merical evidence 11, 131-1611 suggests that plasmoid re- 
connection is "fast", i.e., independent of S. 

In this Letter, we propose a conceptual model of a re- 
sistively reconnecting incompressible plasmoid chain and 
infer the following basic properties: (i) the global recon- 
nection rate is independent of S; (ii) the chain is de- 
scribed by a power law (inverse square) distribution of 
plasmoid fluxes and sizes; (iii) there is a finite probability 
of abnormally large plasmoids that can disrupt the chain. 
We also provide a criterion for a transition from the fast 
resistive MHD regime to a faster collisionless regime. 



Plasmoid chain: general physical picture and key as- 
sumptions. We envision the plasmoid-dominated recon- 
nection layer in a statistical steady state as a chain of 
plasmoids of various sizes separated by small current 
sheets. There is an underlying velocity gradient along 
the chain: all plasmoids are moving outwards as the re- 
connecting flux is coming in. The plasmoids differ greatly 
in size and hence should not be regarded as equal mem- 
bers of the chain — the chain has a hierarchical struc- 
ture. Consider the global current layer of half-length L. 
We may call the largest plasmoids in it the secondary 
ones. Any pair of two adjacent secondary plasmoids 
transforms the region between them into a secondary 
reconnection layer of length 2L^ <C 2L with end-to- 
end longitudinal velocity difference of ~ 2Va- Then, the 
characteristic ejection time for tertiary plasmoids cre- 
ated by the plasmoid instability Q in this secondary 
layer is ~ = /Va, much shorter than the global 
Alfven time ta = L/Va- Thus, the tertiary plasmoids are 
ejected from the secondary layer relatively rapidly and 
their typical size remains smaller than that of the sec- 
ondary plasmoids. Upon ejection, the tertiary plasmoids 
coalesce with the larger secondary ones, contributing to 
their growth. Thus, the secondary layer itself becomes 
a plasmoid-dominated reconnection layer. A similar ar- 
gument can be applied to the tertiary layers and so on. 
We thus get a self-similar hierarchy of plasmoids and in- 
terplasmoid current layers that are themselves complex 
chains of next-generation plasmoids [l4 , li|, 12 1 ■ 



The smallest elementary structure in the chain is the 
"critical layer" — an SP layer marginally stable to the 
plasmoid instability (S = S c ~ 10 ). Its key parameters 
— the length L c = S c r//VA, the thickness S c = L c /y/S^, 
and the reconnection rate cE c = BqVa / a/S'c — depend 
only on n and Va but not on the system size L [lij ]. 
We expect that the smallest current layers found in the 
system are never much shorter or much longer than the 
critical length L c and their reconnection rate always hov- 
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crs around E c . Indeed, because of the underlying ve- 
locity gradient, the secondary current layer between any 
two adjacent plasmoids is continuously stretched. When 
its length exceeds L c , this layer becomes slightly super- 
critical to the plasmoid instability [|| and produces a new 
plasmoid flanked by two new X-points. This new plas- 
moid grows rapidly and reaches the critical size at which 
the two X-points on its sides undergo an X-point col- 
lapse [1,0], turning promptly into two new current sheets. 
This process repeats continuously. We assume that this 
effective splitting of a super-critical inter-plasmoid layer 
into two occurs faster than the typical stretching time of 
that layer; then the lengths of the resulting new layers 
are always somewhat smaller than L c . Thus, the plas- 
moid chain contains N ~ L/L c ~ S/S c plasmoids (of 
all gene rations) separated by nearly critical current lay- 




ers [23[. The hierarchy of plasmoid-dominated layers is 
truncated at the critical layer. 

We can now formalize the above physical picture in 
terms of the following assumptions (or conjectures), ex- 
pected to hold on average: 

(I) The X-point collapses and layer-splitting instabili- 
ties are sufficiently fast (Alfvenic or super- Alfvenic), so 
on average any two immediately adjacent plasmoids are 
separated by a critical layer. 

(II) The upstream reconnecting field in each interplas- 
moid layer is equal to the global reconnecting field B . 
This implies that the outflows into all plasmoids are 
Alfvenic, with velocity Va- This is easy to show via a 
theoretical analysis of the acceleration of plasmoids along 
the layer by the magnetic tension of the associated open 
flux [I9( , but here will be assumed without proof. It then 
follows that the mean flow in the global layer (and in each 
sub- layer) is roughly Hubble- like, v y ~ Vau/L (x and y 
are the directions across and along the layer). 

(III) Plasmoids do not saturate before they are ejected 
from the current layers in which they are embedded into 
larger plasmoids flanking these layers. 

Assumptions I-II appear to be supported by numerical 
evidence but need to be checked systematically. Assump- 
tion III will be verified a posteriori in our theory. 

Effective reconnection rate. A key question in any re- 
connection study is that of the reconnection rate. Be- 
cause of the inherently non-steady nature of plasmoid 
reconnection, we are interested in the effective time- 
averaged (global) rate of transfer of magnetic flux from 
the upstream to the downstream region. To find it, con- 
sider a simple chain of secondary plasmoids and current 
sheets. The fact that the secondary current layers in 
the hierarchical picture are not simple SP sheets but, in- 
stead, complex plasmoid-dominated layers is not impor- 
tant here; we will just view them as effective reconnection 
regions with reconnection rate E^ 2 > . Our goal is to relate 
the effective global reconnection rate E e R to E^> . 

Note that a fully closed magnetic island carries no net 
reconnected flux (B x ). The only reconnected flux actu- 



FIG. 1: Contour plot of 9 from an MHD simulation at S = 
3 ■ 10 5 (detail), illustrating the open magnetic flux. 

ally carried out is the open flux threading the midplane 
x = between the secondary current sheets and plas- 
moids and ejected from the layer with them. This flux 
accumulates because the plasmoid chain is not quite uni- 
form in y: the reconnection rates at the two X-points on 
both sides of any given plasmoid slightly differ (because 
Bq decreases outward along y), resulting in the growth 
of open flux on the faster-reconnecting (smaller- \y\) side 
of the plasmoid. Thus, strictly speaking, the plasmoid 
chain is not a simple sequence of closed magnetic islands 
and current sheets. For example, while the left Y-point 
of a secondary layer located to the right of the global 
center separates it from a plasmoid to its left, the right 
Y-point separates it from a region of open flux followed 
by the plasmoid to its right (see Fig. [I]). 

We can now relate the open flux to the overall recon- 
nection rate. The total reconnected flux is the sum of all 
open-flux parcels between the center of the global layer 
(y = 0) and its right end (y = +L). Reconnection in any 
given secondary layer adds to the open flux on its right 
and subtracts from the open flux enveloping the plasmoid 
on its left. Upon summing up over all the y > layers, 
the contributions from all except for the very first one 
(closest to the global center) cancel. The total open flux 
is thus equal to the flux reconnected via that central- 
most layer. Since in a statistical steady state each plas- 
moid and each open-flux parcel eventually get ejected, 
the net effective reconnection rate equals that of a single 
secondary layer: E g ~ E^ 2 \ Note that this argument is 
very general; it is purely kinematic and does not rely on 
any dynamical considerations (e.g., on our assumptions I- 
III). Applying the same reasoning to all levels n in the 
self-similar plasmoid hierarchy, down to the critical layer 
at the bottom, we see that the effective reconnection rate 
is scale invariant (the same for all levels) and is equal to: 



J eff 



cE c ~ fi-Va v A Bo ~ 1Q- 2 VaB q . (1) 



This means that truly slow reconnection cannot occur, 
even in collisional plasmas. And indeed, we arc not aware 
of any numerical or observational evidence for E c s = 
cE e s/BoVA much below I0~ 2 . 

Plasmoid growth and flux distribution. We are inter- 
ested in a more detailed statistical description of the hi- 
erarchical plasmoid chain, beyond just the effective re- 
connection rate. Namely, we want to find the plasmoid 
spectrum, i.e., the plasmoid flux and size distributions 
(c.f. jl^|; see [l6| for a numerical study of the statis- 
tics of current sheets). To do this, we analyze the inter- 
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play between plasmoid growth, ejection, and coalescence. 
We first focus on "normal" plasmoids — those born not 
very close to the centers of their host layers (previous- 
generation layers in which they are embedded). The im- 
portance of this distinction will become clear later. 

Consider a normal plasmoid of class n inside a host 
layer of length 2L( n_1 ) connecting two plasmoids of 
class n — 1. The plasmoid lives until it is ejected from 
the host layer and merges (coalesces) with a larger plas- 

(n) 

moid of a previous class. The ejection lifetime is t^ ~ 

Z/™ -1 )/^™^, where is the typical plasmoid veloc- 

ity within the host layer. Since, by corollary to as- 
~ Va for all plasmoid classes fljj ]. 

r ]r ^ = -k*'™ -1 ) /Va- During this time, the plas- 
moid grows via reconnection in the two nth-level layers 
on its sides. By assumption III (to be verified later), the 
plasmoid is ejected before saturating, so its flux grows as 



sumption II, V^} 



An) 



tflfW/dt = cE, 



(n+1) 
iff 



cE e s, and eventually reaches 



~ cE ^ +1) t W „ cEcS L^/Va ~ E^BoL^. 

We can now find the cumulative plasmoid-flux distri- 
bution function 7V(\I/) — the total number (in the global 
layer of length 2L) of plasmoids with fluxes larger than \P. 
On the one hand, the relationship between the typi- 
cal flux ty( n ' of normal plasmoids and their host- layer 
length ZA™ -1 ) is given by Eq. ([2]). On the other hand, 
j g comparable to the typical separation between 
plasmoids larger than (i.e., of all previous classes): 
2L(™- 1 ) - Aj/("- x ) - 2L/N{^). Dropping the class 
index n, we then get *iV(*) — cE cS L/Va = E qS B L, 
and hence the flux distribution density is 



/(*) = - dN/dV ~ E cff B L^- 



(3) 



This is a testable quantitative result. Another way to de- 
rive it is to consider a plasmoid-dominated layer and ran- 
domly pick a plasmoid with some flux ^o- Its expected 
age is Tpast ~ ^o/ cE e s and its future life expectancy (the 
time before it is ejected into a larger plasmoid) is Tf utU re ~ 
Ay($ > ^ )/Va ~ L/V A N(^ ). Since the plasmoid was 
chosen randomly, we expect Tf utU re ~ Tpast poT ]- which 
again yields * -W(*o) ~ cE cfi L/V A - 

Growth and distribution of plasmoid sizes. When a 
small amount of flux <5\ff(") = Bq8x is reconnected and 
added to a growing plasmoid of class n, its area is in- 
creased by SA^I ~ Ay^Sx = Ay^5^/B , where 
A?/™) is the separation between plasmoids of class n. On 



the normal-plasmoid life timescale t 



O) 



< 



>-i) 



is not changed strongly by the Hubble flow, Ay^ ~ 
const. Then, the plasmoid area grows just as A' n ' ~ 
Ay("'vl/( ri ) /Bq. As long as the plasmoid x- width w x re- 
mains smaller than its y-extent Wy = Ay( n ) — 2L^ n \ the 
latter stays roughly constant and comparable to Ay^ n \ 
Then, the growth of the plasmoid area translates directly 



into the growth of its a;- width: w x 1 ^ ~ / w : 
*( n V-Bo- Using Eq. ([2j), we find 



00 



,00 



p (rt+l) An) 
C£/ cff l ej 

Bn 



R T/( n ) 



E eS L^ 



(4) 

Interestingly, this allows us to make a comparison with 
the Shibata-Tanuma [171 estimate for the reconnection 



rate cE 



(ST) 
eff 



V out B w max /L, based on a SP mass con- 



servation argument and the assumption that the effective 
outflow channel width is the plasmoid chain width iu roax - 
Applying their estimate to any nth level in the hierar- 
chy and taking w ma 



cE 



(n) 



w x n \ given by Eq. Q, we get: 



V™ Bow^/L^- 1 ) ~ cE^ +1) , independent 
of V^ ut or of n, which coincides with our result Eq. (TTJ). 

The size distribution of normal plasmoids now follows 
from Eq. (|3J) and the width-flux relation w x (i&) ~ ^/B : 



/K) = -dN(w x )/dw x ~ fieifi^ 2 



(5) 



LacA; o/ saturation. These results allow us to verify 
our assumption III that plasmoids do not saturate. Satu- 
ration would occur if w x at any level became comparable 
with the typical separation Ay(^) ~ L/N(^) between 
plasmoids of this size or larger. Using w x ~ ^ / Bq and 
Eq. ©, we have w x /Ay(V) ~ ^N{^)/B Q L ~ E eS < 1, 
Thus, the nonlinear saturation never becomes an issue. 

It is worth noting the crucial role of plasmoid coa- 
lescence (ejection into larger plasmoids) in mitigating 
nonlinear saturation. If we only had a simple chain of 
N ~ L/L c ~ S/S c plasmoids moving along the global 
layer, then plasmoid growth would saturate at w sa t ~ L c , 
corresponding to a plasmoid flux of ^ sat ~ BqL c . The 

1 /2 

time for this to happen, t sat ~ * 8at /aB c ~ (L c /Va)S c ' , 
would be shorter than the global ejection time ta = L/Va 

3 /2 

if S > S c ~ 10 . Thus, without coalescence, plas- 
moids would quickly saturate, stifling reconnection for 

3/2 

S > S c ■ Coalescence prevents the accumulation of sat- 
urated plasmoids and allows larger plasmoids to grow in 
size by eating up smaller ones; it is thus essential for the 
fast plasmoid-dominated reconnection. 

Monster plasmoids. The above picture is modified by 
the presence, in addition to the "normal" plasmoids, of 
relatively rare "anomalous" plasmoids born near the cen- 
ters of their host layers at any level in the hierarchy. 
The most important of them are the lowest-generation 
anomalous plasmoids born near the center of the global 
layer of length 2L. Because the Hubble-flow ejection 
time for an anomalous plasmoid born at y = yo <C L 
is i C j = f dy/vy ~ ta log(£/yo)j this plasmoid lives 
longer than normal plasmoids and hence grows larger. 
Its final flux is enhanced just by a logarithmic factor: 
\& ~ cE e gTA log(L/yo); however, the enhancement of its 
area is much greater. Indeed, a given plasmoid grows 
in mass and area by sucking in all the plasma (includ- 
ing smaller plasmoids) within its domain of influence, 
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which extends up to the mid-point between it and the 
next plasmoid of similar size. Thus, the area growth rate 
is proportional to the interplasmoid separation Ay(t): 
dA/dt ~ [Ay(t)/B }dV/dt ~ Ay(t) E eS V A - But, on 
their long ejection timescale, the separation between 
two anomalous plasmoids grows exponentially, Ay(t) ~ 
Ay(0) ex.p(t/r A ). Therefore, A(t) also grows exponen- 
tially and ultimately reaches A max = f " dt (dA/dt) ~ 
E e g L 2 Ay(0)/i/o, larger by a factor L/ya ^> 1 than for 
normal plasmoids. Taking the typical smallest initial po- 
sition j/o to be ~ Ay(0), we get A max ~ E e gL 2 . Next, it 
is easy to show that for S > Sc^ 4 ~ 10 5 , this growth is so 
rapid that the plasmoid a;-width w x catches up with its 
initial y-extent w y ~ Ay(0) before t C j is reached. Then, 
subsequently, w x and w y grow in unison and eventu- 



ally reach w B 



A 



1/2 



Ell 2 L ~ 0.1 L. This pre- 



diction of occasional large, macroscopic "monster" plas- 
moids should have important implications for observa- 
tions (e.g., large abrupt events in solar flares [21( and 
sawtooth crashes 0) and simulations. 

Transition to collisionless reconnection. Even if the 
global reconnection layer is in the resistive MHD regime, 
8sp{L) > min{dj, p s }, where di is the ion collisionless 
skin-depth and p s is the ion sound Larmor radius, this 
may not be so for the smaller layers in the plasmoid hier- 
archy. Namely, if S c < minjdi, p s }, then resistive MHD 
breaks down at some level in the hierarchy and the cor- 
responding current sheets transition to fast collisionless 
reconnection with cEnaii — OABqVa ~ 10 cE r .. The hier- 
archy terminates at this point (c.f. [12, 3, 22|). Most of 
our results should still hold, with E c g ~ ^Hali- Interest- 
ingly, the resulting range, 0.01 < E e g < 0.1, covers most 
of the rates inferred observationally or numerically. 

One can evaluate the ratio S c /di as 



s c /di ~ sl' 2 v/v A di ~ sl' 2 v e /n e 



(6) 



where f2 e and v e are the electron cyclotron and colli- 
sion frequencies (c.f. (l2|). Thus, the plasmoid hierar- 
chy stays collisional all the way down to the critical layer 
only in very collisional plasmas. Using Spitzcr resistivity, 
S c /di ~ (5 c 1/2 /12tt) (c/V a ) (me/mA 1 / 2 lnAiV" 1 , where 
Nr> = (47r/3) n e A|j is the number of electrons in a Debye 
sphere. Since Nu 3> 1 for a medium to be considered a 
normal plasma, we expect S c <C di and so the smallest 
layers in the plasmoid chain are unavoidably collisionless. 
In the strong-guidc-ficld case, di is replaced by p s and the 
above expressions are only changed by factors of ~ ft 1 ' 2 , 
which does not significantly affect our conclusions. 

Conclusions. We have proposed a simple model of re- 
connection in plasmoid dominated current layers which 
yields a Lundquist-number independent effective recon- 

— 1/2 

nection rate E e g ~ S c ~ 0.01 and a self-similar dis- 
tribution of plasmoid sizes and fluxes (an inverse-square 
law, a testable prediction). In addition, we have argued 
that the plasmoid instability induces a multi-level plas- 



moid hierarchy that almost always reaches the kinetic 
scales, implying that pure resistive MHD reconnection 
occurs only in the most collisional plasmas. Our pre- 
diction of the occurrence of monster plasmoids, of width 
w x ~ 0.1L, offers a possible interpretation of observable 
events in solar flares and sawtooth crashes. 

While a detailed nonlinear theory of plasmoid recon- 
nection remains a challenge, as do fully resolved simula- 
tions of such a process, we hope that the simple model 
presented above might provide a useful conceptual frame- 
work for high-Lundquist-number reconnection in the way 
the SP model has done for moderate Lundquist numbers. 
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